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ABSTRACT 

The nova-Hke cataclysmic binary AE Aqr, which is currently understood to be a former supersoft X-ray 
binary and current magnetic propeller, was observed for over two binary orbits (78 ks) in 2005 August with the 
High-Energy Transmission Grating (HETG) onboard the Chandra X-ray Observatory. The long, uninterrupted 
Chandra observation provides a wealth of details concerning the X-ray emission of AE Aqr, many of which are 
new and unique to the HETG. First, the X-ray spectrum is that of an optically thin multi-temperature thermal 
plasma; the X-ray emission lines are broad, with widths that increase with the line energy, from ct « 1 eV 
(510 km s"') for O VIII to a k, 5.5 eV (820 km s"') for Si XIV; the X-ray spectrum is reasonably well fit 
by a plasma model with a Gaussian emission measure distribution that peaks at logr(K) = 7.16, has a width 
a = 0.48, an Fe abundance equal to 0.44 times solar, and other metal (primarily Ne, Mg, and Si) abundances 
equal to 0.76 times solar; and for a distance d = 100 pc, the total emission measure EM = 8.0 x 10^^ cm~^ and 
the 0.5-10 keV luminosity Lx = 1-1 x 10^' erg s"'. Second, based on the f/{i + r) flux ratios of the forbidden 
(/), intercombination (/), and recombination (r) lines of the He a triplets of N VI, O VII, and Ne IX measured by 
Itoh et al. in the XMM-Newton Reflection Grating Spectrometer spectrum and those of O VII, Ne IX, Mg XI, and 
Si XIII in the Chandra HETG spectrum, either the electron density of the plasma increases with temperature 
by over three orders of magnitude, from We ~ 6 x 10'" cm"^ for N VI [log r(K) « 6] to We ~ 1 x 10''* cm"^ for 
Si XIII [log r(K) « 7], and/or the plasma is significantly affected by photoexcitation. Third, the radial velocity 
of the X-ray emission lines varies on the white dwarf spin phase, with two oscillations per spin cycle and 
an amplitude K « 160 km s"'. These results appear to be inconsistent with the recent models of Itoh et al., 
Ikhsanov, and Venter & Meintjes of an extended, low-density source of X-rays in AE Aqr, but instead support 
earlier models in which the dominant source of X-rays is of high density and/or in close proximity to the white 
dwarf. 

Subject headings: binaries: close - stars: individual (AE Aquarii) - novae, cataclysmic variables - X-rays: 
binaries 



I. INTRODUCTION 

AE Aqr is a bright (V « 11) nova-like cataclysmic binary 
consisting of a magnetic white dwarf primary and a K4-5 V 
secondary with a long 9.88 hr orbital period and the short- 
est kn own white dwarf spin period P = 33.08 s (iPattersonI 
119791) . Although originally classi fied and i nterpreted as a 
disk-accreting DQ Her star (Patte rsonI 11994 1). AE Aqr dis- 
plays a number of unusual characteristics that are not natu- 
rally explained by this model. First, violent flaring activity 
is observed in the radio, optical, ultraviolet (UV), X-ray, and 
TeV 7-rays. Second, the Balmer emission lines are single- 
peaked and produce Doppler tomograms that are not consis- 
tent with those of an accretion disk. Third, the white dwarf is 
spinning down at a rate P = 5.64 x lO"'** s s"' dde Jager et al.l 
[1994) . Although this corresponds to the small rate of change 
of 1 .78 ns yr"' , AE Aqr's spin-down is typically characterized 
as "rapid" because the characteristic time P/P « 2 x 10^ yr is 
short compared to the lifetime of the binary and because the 
spin-down luminosity Lsd = -I^^ w 1 x lO^'* erg s"' (where 
/ « 0.2Mwd/?^^ « 2 x 10^" g cm^ is the moment of inertia for a 
white dwarf of mass M^d = 0.8 Mq and radius R^d = 7.0 x 10^ 
cm, n = 2tt/P, and Ct = -2ttP/P^) exceeds the secondary's 
thermonuclear luminosity by an order of magnitude and the 
accretion luminosity by two orders of magnitude. 

Because of its unique properties and variable emis- 
sion across the electromagnetic spectrum, AE Aqr has 
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been the subject of numerous studies, including an inten- 
sive multiwayelengt h observing campaign in 1993 October 
llCasares et al.lll996i and the serie s of papers in ASP Conf. 
Ser. 85 (iBucklev & W arner '1995)1. Based on these stud- 
ies, AE Aqr is now widely believed to be a former super- 
soft X-ray binary (Schenker et al. 2002) and current mag- 
netic propeller (Wynn, King, & Ho rne 1997), with most of 
the mass lost by the secondary being flung out of the bi- 
nary by the magnetic field of the rapidly rotating white 
dwarf. These models explain many of AE Aqr's unique 
characteristics, including the fast spin rate and rapid sec- 
ular spin-down rate of the white dwarf, the anomalous 
spectral type of the secondary, the anomalous abundances 
dMauche. Lee. & K allman 1997,). the absen ce of signatures of 
an accretion disk (Welsh. Horne. & Gomeii[T9 98). the violent 
flaring activity ([Pearson. Horne. & Skidm ore 2003), a nd the 
orig in of the ra dio and TeV 7-ray emission (Kuiiper s et al.l 
1997; Meintjes & Venter 2003). 

To build on this observational and theoretical work, while 
taking advantage of a number of improvements in observing 
capabilities, during 2005 August 28-September 2, a group of 
professional and amateur astronomers conducted a campaign 
of multiwavelength (radio, optical, UV, X-ray, and TeV 
7-ray) observations of AE Aqr. Attention is restricted 
here to the results of the X-ray observations, obtained with 
the High-Energy Transmission Grating (HETG) and the 
Advanced CCD Imaging Spectromete r (ACIS) dete ctor 
onboard the Chandra X-ray Observatory. iMauchd ( 120061) has 



2 



Mauche 



previously provided an analysis and discussion of the timing 
properties of these data, showing that: (1) as in the optical 
and UV, the X-ray spin pulse follows the motion of the white 
dwarf around the binary center of mass and (2) during the 
decade 1995-2005, the white dwarf sp un down at a r a te tha t 
is slightly faster than predicted by the |d e Ja ger et al j (119941) 
spin ephemeris. Here, we present a more complete analysis 
and discussion of the Chandra data, providing results that in 
many ways reproduce the results of th e previous Einstein , 
ROSAT, ASCA, and XMM-Newton (Patterson et al.' '19801 
Reinsch et al. 1995; Clavton & Osborne 1995; Osborne et 
T995i; [EracleousI 11 999; Choi, Dotan i. & Agrawall \\99^ 
Itoh et alj I2006t IChoi & Dotani 2006.) and the subsequent 
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Suz.aku dTerada et alJ 120081) observations of AE Aqr, but 
also that are new and unique to the HETG; namely, the 
detailed nature of the time-average X-ray spectrum, the 
plasma densities implied by the He a triplet flux ratios, and 
the widths and radial velocities of the X-ray emission lines. 
As we will see, thes e results appear t o be inconsistent with 
the recent models of lltoh et alj (120061) , llkhsanovl (120061) , and 
I Venter & MeintiesI (l2007h of an extended, low-density source 
of X-rays in AE Aqr, but instead support earlier models in 
which the dominant source of X-rays is of high density and/or 
in close proximity to the white dwarf. 

The plan of this paper is as follows. In §2 we discuss the 
observations and the analysis of the X-ray light curve (§2.1), 
spin-phase light curve (§2.2), spectrum (§2.3), and radial ve- 
locities (§2.4). In §3 we provide a summary of our results. 
In §4 we discuss the results and explore white dwarf (§4.1), 
accretion column (§4.2), and magnetosphere (§4.3) models of 
the source of the X-ray emission in AE Aqr In §5 we draw 
conclusions, discuss the bombardment model of the accretion 
flow of AE Aqr, and close with a few comments regarding fu- 
ture observations. The casual reader may wish to skip §2.1- 
2.2, which are included for completeness, and concentrate on 
§2.3-2.4, which contain the important observational and anal- 
ysis aspects of this work. 

2. OBSERVATIONS AND ANALYSIS 

AE Aqr was observed by Chandra beginning on 2005 Au- 
gust 30 at 06:37 UT for 78 ks (ObsID 5431). The level 2 
data files used for this analysis were produced by the stan- 
dard pipeline processing software ASCDS version 7.6.7. 1 and 
CALDB version 3.2.1 and were processed with the Chan- 
dra Interactive Analysis of Observations (CIAO') version 3.4 
software tools to convert the event times in the evt2 file from 
Terrestrial Time (TT) to Barycentric Dynamical Time (TDB) 
and to make the grating response matrix files (RMFs) and the 
auxiliary response files (ARFs) needed for quantitative spec- 
troscopic analysis. The subsequent analysis was performed 
in the following manner using custom Interactive Data Lan- 
guage (IDE) software. First, using the region masks in the 
pha2 file, the source and background events for the ± first- 
order Medium Energy Grating (MEG) and High Energy Grat- 
ing (HEG) spectra were collected from the evt2 file. Sec- 
ond, after careful investigation, ±0.0030 A (±0.0015 A) was 
added to the ± first-order MEG (HEG) wavelengths, respec- 
tively, to account for an apparent shift (by 0.27 ACIS pix- 
els) in the position of the zero-order image. Third, to ac- 
count for the spin pulse delay measured in the optical, UV, and 
X-ray waveban ds dde Japer et al.lll994tlEracleous et alJI 19941 ; 
iMauchd [20061) produced by the motion of the white dwarf 
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Fig. 1. — MEG plus HEG count rate light curve of AE Aqr {black his- 
togram) and 1(7 en'or vector (lower gray histogram). Bin width At = 300 
s. 

around the binary center of mass, -2 cos 27r(/)orb s was added to 
the event times f , where the white dwarf orbit phase 27r(/iorb = 
f^orb(f-7()), where rJorb = 27r/Porb and P^a, = 0.411655610 d 
and r()(BJD) = 244517 2.2784 are the orbit ephemeris con- 
stants from Table 4 of Ide Jager et al.l (jT994). Fourth, to ac- 
count for the Doppler shifts produced by Chandra's (mostly. 
Earth's) motion relative to the solar system barycenter, the 
event wavelengths A were multiplied by a factor [1 + vios/c], 
where vios is the (time-dependent) line-of-sight velocity be- 
tween the spacecraft and the source, determined using an 
IDE code kindly supplied by R. Hoogerwerf (and checked 
against the line-of-sight velocities derived from the barycen- 
tric time corrections supplied by the CIAO tool axbary); 
during the observation, vios varied from -11.9 km s"' at the 
beginning of the observation, rose to -11.3 km s"', and then 
fell to -11.5 km s"' at the end of the observation. Fifth, 
a filter was applied to restrict attention to events from two 
orbital cycles 0orb = 20503.9-20505.9, resulting in an effec- 
tive exposure of 71 ks. Sixth, white dwarf spin pulse phases 
were calculate d using the updated cubic spin ephemeris of 
IMauchd (120061) derived from the recent A5CA (1995 October) 
and XMM-Newton (2001 November) and the current Chan- 
dra observation of AE Aqr: liKp^^xj^ = fio(f - Tmax) + 
7max)^+ 5fj(f- Tmax)^ whcrc Qq = 27r/f 33 and n = -InPii/Pl 
and P33 = 0.00038283263840 d, P33 = 5.642 x lO"'"* d d"', 
and rinax(BJD) = 2445 172.000042 are the spin ephemeris con- 
stants from Table 4 of Ide Jager et all (119941) , and f2 = -1 .48 x 
10~" d"^. The HETG event data are then fully characterized 
by the event time t, white dwarf orbit phase ^orb, white dwarf 
spin phase 0spin, and wavelength A. 

2.1. Light curve 

Figure 1 shows the background-subtracted count rate light 
curve derived from the HETG event data. As has been well 
established by previous X-ray observations, the X-ray light 
curve of AE Aqr is dominated by flares, although this is by 
far the longest uninterrupted observation of AE Aqr and hence 
the clearest view of its X-ray light curve. During the Chandra 
observation, the flares last between a few hundred and a few 
thousand seconds, producing increases of up to 3-5 times the 
baseline count rate of 0. 1 counts s"' . 

To constrain the cause and nature of the flares of AE Aqr, 
it is of interest to determine if the count rate variations shown 
in Figure 1 are accompanied by, or perhaps are even due to. 
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Fig. 2. — Soft (1 1-26 A) over hard (1-8 A) versus total (1-26 A) count rate 
for AE Aqr {filled circles with error bars) and best-fitting constant function 
{dotted line) for the first three data points. 

dramatic variations in the X-ray spectrum. Previous investi- 
gations have indicated that this is not the case. For a flare ob- 
served by ASCA, Ch oi. Do t ani, & Agrawal (1999) found no 
significant difference between the quiescent and flare X-ray 
spectra, although a "hint" of spectral hardening was recog- 
nized during th e flare. For a flare observed by XMM-Newton, 
IChoi & Dotanil (12006 ^ found that the X-ray spectrum at the 
beginning of the flare was similar to that in quiescence, but 
that the spectrum became harder as the flare advanced. Our 
ability to investigate spectral variations during the individ- 
ual flares observed by Chandra is limited by the relatively 
low HETG count rate and the relatively fast timescale of the 
flares, although it is possible to investigate spectral variations 
for the ensemble of the flares. To accomplish this, the light 
curve shown in Figure 1 was divided into five count rate 
ranges: h < 0.14 counts s"', 0.14 < /a < 0.18 counts s"', 
0.18 < /3 < 0.23 counts s"', 0.23 <h< 0.29 counts s"', 
and /5 > 0.29 counts s~', where the count rate cuts were set 
to produce a roughly equal number of counts per count rate 
range, and the source and background counts were collected 
in three wavebands: hard (1-8 A), medium (8-1 1 A), and soft 
(11-26 A), where the wavelength cuts were set to produce a 
roughly equal number of counts per wavelength interval. The 
background-subtracted soft S over hard H versus the total (1- 
26 A) count rates are plotted in Figure 2. The dotted line in 
that figure is the best-fitting constant function 5/// = 0.953 for 
the first three data points (count rate ranges /1-/3). As shown 
by the figure, the softness ratio of the next higher count rate 
range I4 is consistent with the lower ranges, while that of the 
highest count rate range A is significant l y (3.6 g) less. Con- 
sistent with the result of lChoi & Dotanil (l2006h . we find that 
the X-ray spectrum of the flares of AE Aqr is harder than that 
in quiescence only at the peak of the flares. 

2.2. Spin-phase light curve 

Figure 3a shows the background-subtracted spin-phase 
count rate light curve derived from the HETG event data. It 
is well fit (x^ per degree of freedom = = 7.26/7 = 1.04) 
by the cosine function A + B cos27r((/)spin - (/)o), with mean 
count rate A = 0.191 ±0.002 counts s"', semi-amplitude B = 
0.031 ±0.002 counts s"' (hence, relative X-ray pulse semi- 
amplitude B/A = 16% ±1%) a nd, consistent with the updated 
spin ephemeris ( lMauchell2006l) . phase offset 4>q = 0.00 ± 0.01 
(throughout the paper, errors are \a or 68% confidence for 
one free parameter). 

The above result appUes to the observation-average spin- 
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Fig. 3. — {a) Spin-phase count rate light curve of AE Aqr {filled circles with 
error bars), best-fitting cosine function {solid curve), and mean count rate A 
{dotted line), (b) Soft (11-26 A) over hard (1-8 A) spin-phase count rate 
light curve of AE Aqr {filled circles with error bars) and best-fitting constant 
function {solid line). 
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Fig. 4. — Semi-amplitude B versus mean count rate A for the intensity- 
resolved spin-phase count rate light curves of AE Aqr {filled circles with error 
bars), best-fitting constant fit to B/A {dotted line) and best-fitting linear fits 
to B/A versus A {dashed curve) and B versus A {solid line). 

phase light curve, although it is of interest to determine if 
the X-ray pulse semi-amplitude B and/or the relative semi- 
amplitude B/A varies with the mean count rate A. As for the 
investigation of the softness ratio variations, this is best de- 
termined as a function of time, using a time resolution suf- 
ficient to resolve the flares, but this is not possible with the 
relatively low count rate HETG data. Instead, background- 
subtracted spin-phase count rate light curves were derived 
from the HETG event data for each of the five count rate 
ranges defined above, and were then fit with the cosine func- 
tion A +Z? cos27r0spin- The resulting values of the X-ray pulse 
semi-amplitude B are plotted versus the mean count rate A 
in Figure 4. The figure demonstrates that B increases linearly 
with A {B/A is constant with A) in the middle of the count rate 
range, but that it saturates at both the low and high ends of the 
range; over the full count rate range, B/A is not well fit (x^ = 
5.82/4= 1 .46) by a constant B/A = 0. 162 ± 0.012 (dotted line 
in Fig. 4). Instead, B/A versus A and B versus A are well fit 
by a linear relation a + bA with, respectively, a = 0.22 ±0.03 
and b = -0.25 ±0.15, with xl = 2.77/3 = 0.92 (dashed curve 
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Fig. 5. — MEG plus HEG count spectrum of AE Aqr (black histogram) and la error vector (lower gray histogram). Bin width AA = 
used to construct the composite line profile are labeled. 
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in Fig. 4), and a = 0.011 ±0.006 and b = 0.10 ±0.03, with 
xl = 2.38 /3 = 0.79 (solid line in Fig. 4). 

To determine if the observed spin-phase flux modulation is 
due to photoelectric absorption (or some other type of broad- 
band spectral variability), background-subtracted spin-phase 
count rate light curves were derived from the HETG event 
data for each of three wavebands defined above: hard (1-8 
A), medium (8-11 A), and soft (11-26 A). Figure 3b shows 
the ratio of the resulting soft S over hard H spin-phase count 
rate light curves. The data is well fit (xl = 7.20/9 = 0.80) 
by a constant S/H = 0.916 ±0.019, which strongly constrains 
the cause of the observed spin-phase flux modulation. Specif- 
ically, if the observed flux modulation is caused by photo- 
electric absorption, a variation in the neutral hydrogen col- 
umn density AA^h ~ 3 x 10^' cm~^ is required, whereas the 
essential constancy of the softness ratio light curves requires 
AA^H < 1 X 10^" cm"^; a factor of 30 times lower. 



2.3. Spectrum 

Figure 5 shows the background-subtracted count spectrum 
derived from the HETG event data, using AA = 0.05 A wave- 
length bins as a compromise between spectral resolution and 
sign-to-noise ratio. As is typical of unabsorbed cataclysmic 
variables (CVs) (Mukai et al. 2003; Mauche 2007), the X-ray 
spectrum of AE Aqr is that of a multi-temperature thermal 
plasma, with emission lines of H- and He-like O, Ne, Mg, Si, 
and S and L-shell Fe XVII-Fe XXIV. However, unlike other 
CVs, and, in particular, unlike other magnetic CVs, the H- 
like Fe XXVI and He-like Fe XXV lines and the "neutral" flu- 
orescent Fe K line are not apparent in the HETG spectrum. 
The apparent absence of these features in the HETG spec- 
trum places limits on the maximum temperature of the plasma 
in AE Aqr and the amount of reflection from the surface of 
the white dwarf, although higher-energy instruments, such as 
those onboard Suzaku, are better suited to study this portion 
of the X-ray spectrum (see Terada et al. 2008). 

We conducted a quantitative analysis of the X-ray spectrum 
of AE Aqr in three steps. First, Gaussians were fitted to the 
strongest emission lines of H- and He-like O, Ne, Mg, and Si 
to determine their radial velocities, widths, and fluxes. Sec- 
ond, a global model was fitted to the X-ray spectrum to con- 
strain its absorbing column density, emission measure distri- 
bution, and elemental abundances. Third, using the emission 
measure distribution and the flux ratios of the He a forbidden 
(/), intercombination (/), and recombination (r) lines, con- 
straints are placed on the density of the plasma. 
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Fig. 6. — Gaussian width cr versus line energy for the Lyman a einission 
lines of H-like O, Ne, Mg, and Si (filled circles with error bars) and best- 
fitting linear function (solid line). 

2.3.1. Line radial velocities, widths, and fluxes 

The radial velocities, widths, and fluxes above the contin- 
uum of the Lyman a emission lines of O VIII, Ne X, Mg XII, 
and Si XIV and the He a triplets of O VII, Ne IX, Mg XI, and 
Si XIII were determined by fitting the flux in the MEG spec- 
trum in the immediate vicinity of each emission feature with a 
constant plus one (Lyman a) or three (He a) Gaussians, em- 
ploying the ARE and RMF files to account for the effective 
area of the spectrometer and its AA = 0.023 A [690 km s"' at 
10 A] FWHM spectral resolution. For each emission feature, 
the radial velocity was determined relative to the laboratory 
wavelengths from the Interactive Guide for ATOMDB ver- 
sion 1.3. More specifically, the assumed wavelengths for the 
Lyman a lines are the mean of the wavelengths of the dou- 
blets weighted by their relative emissivities (2: 1), whereas the 
wavelengths for the He a intercombination lines are the un- 
weighted means of their component x and y lines. In the fits 
to the He a triplets, the radial velocities and widths deter- 
mined from the fits to the corresponding Lyman a lines were 
assumed, so that, in all cases, the fits had four free param- 
eters. For these fits, u nbinned data were employed and the 
C statistic jCashl [T979h was used to determine the value of 
and error on the fit parameters, which are listed in Table 1 . 
As demonstrated in Figure 6, we find that the widths of the 
Lyman a emission lines increase with the line energy, from 
o- = 1.1 ±0.1 eV for O VIII to cr = 5.5 ±0.9 eV for Si XIV. 
For comparison, during two flares of AE Aqr observed with 
the XMM - Newto n Reflection Grating Spectrometer (RGS), 
lltoh et alJ (12006b found cr w 1 .2 and 2 eV for the Lyman a 

2 Available at |http://cxc.harvard.edu/atomdbAVebGUIDE/l 



Chandra HETG Spectrum of AE Aqr 



5 



TABLE 1 

Line Fit Parameters, Flux Ratios, and Inferred Electron Densities. 





velocity 


cr 




Flux (10^ photons cm - s ') 








log He 


Element 


(kms-') 


(eV) 


Lyman a 


/ 


(' 


r 


G = (/+OA 


R'=//(! + r) 


(cm-3) 


O 


-75 ± 51 


1.12 ±0.09 


2.91 ±0.27 


0.35 ±0.17 


0.84 ±0.24 


1.12 ±0.26 


1.05 ±0.35 


0.178 ±0.090 


11 -36^32 

12 30+^-^' 

-Q.36 

1 3 nA-KiXS 

-fl.23 

14 M+O-^' 


Ne 


+45 ± 53 


2.02±0.17 


0.999 ± 0.073 


0.206 ±0.067 


0.200 ± 0.064 


0.622 ± 0.089 


0.65 ±0.18 


0.251 ±0.088 


Mg 


+49 ±104 


3.85 ±0.60 


0.188±0.019 


0.070 ±0.021 


0.040 ± 0.020 


0.193 ±0.024 


0.57±0.17 


0.300 ±0.099 


Si 


+95 ±118 


5.50 ±0.90 


0.185±0.019 


0.055 ±0.014 


0.028 ±0.017 


0.205 ±0.020 


0.41±0.11 


0.235 ±0.065 



emission lines of N VII and O VIII, respectively. The trend 
shown in Figure 6 is well fit (x^ = 1.43/2 = 0.71) with a linear 
function a + bE with a = -0.80 ± 0.29 and = 2.9 ± 0.4. Con- 
sistent with the non-zero intercept a, the line widths are not 
constant in velocity units, but increase with the line energy: 
cr = 512 ± 39, 593 ±51, 784 ± 121, and 822 ± 135 km s"' 
for O VIII, Ne X, Mg XII, and Si XIV, respectively. In ad- 
dition to the hne widths, there is some evidence that the ra- 
dial velocities of the Lyman a emission lines increase with 
the line energy, from v = -75 ±51 km s"' for O VIII to 
V = +95 ± 1 18 km s"' for Si XIV. Note, however, that the ve- 
locity difference, Av = 170 ± 128 km s"', differs from zero by 
just 1.3(7, is a small fraction of the widths of the lines, and is 
probably affected by systematic effects. If a common radial 
velocity offset is assumed in the fits of the Lyman a lines, the 
derived velocity v = -1 ± 33 km s~' . 

2.3.2. Global model 

To produce a global model of the X-ray spectrum of AE 
Aqr, ATOMDB IDL version 2.0.0 software^ was used to con- 
struct ATOMDB version 1.3.1"' optically-thin thermal plasma 
X-ray spectral models for the continuum and for the 12 cos- 
mically abundant elements C, N, O, Ne, Mg, Al, Si, S, Ar, 
Ca, Fe, and Ni at 40 temperatures spaced uniformly in log T 
[specifically, logr(K) = 5.0,5.1,5.2, ... ,8.9]. Using custom 
IDL software, these spectral eigenvectors were convolved 
with a Gaussian, to account for the observed widths of the 
emission lines as a function of temperature, and multiplied by 
the grating ARFs and RMFs, to account for the spectrome- 
ter's effective area and spectral resolution. The resulting ± 
first-order MEG and HEG spectral models were then binned 
to 0.05 A and coadded. Finally, the observed MEG plus HEG 
count spectrum (Fig. 5) and whence the spectral models were 
"grouped" to a minimum of 30 counts per bin so that Gaussian 
statistics could be used in the fits. 

Numerous model emission measure {EM) distributions 
were tested against the data: 1, 2, 3, and 4 single-temperature 
models, a cut-off power law (dEM /dlogT cx T" for T < 
Tc), a power law with an exponential cutoff (dEM /dlogT oc 
T" for r < Tc and dEM /dlogT oc T°'ti^p{{%-T)/Tf\ for 
T > Tc), and a Gaussian (dEM /dlogT oc exp[-(logr- 
logr())^ /2cr^]), all with photoelectric ab sorption by a neutral 
column dMorrison & McCammori 1983h and variable elemen - 
tal abundances relative to those of lAnders & Grevessel(ll989t) . 

Among these models, the best (if not a particularly good) 
fit (x^ = 382.7/233 = 1.64) was achieved with a model with 
a Gaussian emission measure distribution with a peak tem- 
perature logro(K) = 7.16 ±0.01, a width a = 0.48 ±0.01, 
an absorbing column density A^h = (1.0 ±0.6) x 10^*^ cm"^, 
an Fe abundance equal to 0.44 ± 0.02 times solar, and the 



^ Available at http://asc.hai'var d.edu/atomdb/features_idl.html| 
"* http://asc.harvard.edii/atomdb/ 



Other metal abundances equal to 0.76 ± 0.03 times solar^ 
The X-ray spectrum of this model is shown superposed on 
the data in Figur e 7. For an assumed distance d = 100 pc 
( lFriediung|ll997h . the total emission measure EM = 8.8 x 
10^^ cm"^ and the unabsorbed 0.5-10 keV luminosity Lx = 
1.1 X 10-" erg s"i; ;/ this luminosity is due to accretion onto 
the white dwarf, the mass accretion rate M = Lx^wd/GM„d w 
7.3 X 10'3gs-i. 

Before leaving this section, it is useful to note that the Gaus- 
sian emission measure distribution derived above extends over 
nearly two orders of magnitude in temperature — from T w 
1 .6 X 10*" K to 1 .3 X 10*^ K at ±2cr — but peaks at a temper- 
ature r « 1.4 X 10^ K or 1.2 keV. Such a low temperature is 
uncharacteristic of CVs, and in particular magnetic CVs: note 
for instance that AE Aq had the lowest continuu m tempera- 
ture o f any magnetic CV observed by ASCA ( Ezu ka & Ishidal 
119991) . The characteristic temperatures of magnetic CVs are 
the shock temperature, T^ = 3fj.miiGM„i/8kR„d « 4.1 x 10^ 
K or 35 keV, applicable for radial free-fall onto the white 
dwarf, and the blackbody temperature, Tbb = (Lx/fA)'/'* w 
19 (//0.25)-i'''* kK (where the radiating area A = AttRIJ and 
our choice for the fiducial value of the fractional emitting area 
/ will be justified in §4.1), applicable for deposition of the 
accretion luminosity in the surface layers of the white dwarf 
(e.g., bombardment or blobby accretion). The only evidence 
for plasma at the shock temperat ure is supplied by the Suzaku 
X-ray spectrum of AE Aqr, but iTerada et al.l ( l2008h argued 
strongly for a nonthermal origin for this emission. The black- 
body temperature, on the o ther hand, is close to the tempera- 
ture of the hot spots that Eracl eous et al.l ('1994') derived from 
the maximum entropy maps of the UV pulse profiles of AE 
Aqr. More on this below. 

2.3.3. Plasma densities 

We now consider in more detail the forbidden (/), in- 
tercombination (i), and recombination (r) line fluxes of 
the He a triplets of O VII, Ne IX, Mg XI, and Si XIII 
der ived in §2.3.1 an d liste d in Table 1. As elucidated 
by [G abriel & JordanI (119691); jBlumenthal. Drake. & Tucked 
( 1972); and Porquet et al. (2001]), these line fluxes can be used 
to constrain the electron temperature T^ via the G= (f+ i)/r 
flux ratio and the electron density rig via the R = f/i flux ra- 
tio of each He-like ion. Because the errors on the observed R 
flux ratios are large, we follow Itoh et al. (2006) and employ 
the R' = f /ii + r) flux ratio as an electron density diagnostic 
in the subsequent discussion. The observed G and R' flux 
ratios and errors are listed in Table 1 . To interpret these re- 
sults, we derived the R' = RG/{1 +R + G) flux ratios from the 
G and R flux ratios tabulated by iPorqu et et al. (2001) for a 

^ This result is driven primarily by, and should be understood to apply 
primarily to, Ne, Mg, and Si. Consistent with the quaUty of the data, no 
attempt was made to further subdivide the element abundances. 
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Fig. 7. — (a) MEG plus HEG spectrum of AE Aqr (black histogram), Icr error vector (lower blue histogram), and the best-fit ATOMDB absorbed variable- 
abundance Gaussian emission measure distribution model (red histogram), (b) Corresponding residuals. 



collisional plasma for N VI, O VII, Ne IX, Mg XI, and Si XIII 
for logre(K) = 6.2, 6.4, 6.6, 6.8, and 7.0, respectively, where 
the temperatures are those of the peaks of the He-like triplet 
emissivities weighted by the Gaussian emission measure dis- 
tribution, determined in the previous section from the global 
fit to the X-ray spectrum. With the exception of Si XIII, * in 
each case the assumed electron temperature is consistent with 
the observed G flux ratio. With these assumptions, the theo- 
retical values of R' are shown in the left panels of Figure 8 as 
a function of \ogn^. In addition to the R' flux ratios of O VII, 
Ne IX, Mg XI, and Si XIII measured from the Chandra HETG 
spectrum (Table 1), we added to Figure 8 the R' flux ratios of 
N VI, O VII, and Ne IX measured by Itoh et al. from the XMM- 
Newton RGS spectrum of AE Aqr. Figure %{b) of this paper 
corrects an error in Figure 5{h) of Itoh et al., which showed 
the R' flux ratio extending from 0.09 to 0.39, whereas the data 
in their Table 4 shows that it should extend only to 0.29. With 
this correction, the RGS- and HETG-derived values of and er- 
rors on the R' flux ratio of O VII are nearly identical. This cor- 
rection, the significantly smaller error range on the R' flux ra- 
tio of Ne IX, and the HETG results for Mg XI and Si XIII indi- 
cates that, contrary to the conclusion of Itoh et al., the electron 
density of the plasma in AE Aqr increases with temperature 
by over three orders of magnitude, from «e ~ 6 x 10'*' cm"-' 
for N VI [log7;(K) w 6.2] to We « 1 x 10'"* cm'^ for Si XIII 
[log7;(K)«7.0]. 

In addition to electron density, the R and hence the R' 
flux ratio is affected by photoexcitation, so we must inves- 
tigate the sensitivity of the R' flux ratios to an external radi- 
ation field. To investigate the conditions under which a low- 
density plasma can masquerade as a high-density plasma, we 
derived the R' flux ratios from the G and R flux ratios tabu- 
lated by Porquet et al. (2001) for a low-density {n^ = 1 cm"-') 
collisional plasma irradiated by a blackbody with tempera- 
ture 7bb and dilution factor W . Under the assumption that 
this flux originates from the surface of the white dwarf, W = 
5 {1 - [l-{Rvji/ r)Y^'^}, where r is the distance from the cen- 
ter of the white dwarf, hence W = ^ on the white dwarf sur- 



In the IPorquet et al] 1200 Ih tabulation, for Si Xin for 



10'- 



G = 0.90, 0.76, 0.67, and 0.56 for T = 5, 7.5, 10, and 15 MK, respectively, 
but then rises to G = 0.70 for 7 = 30 MK; the observed G = 0.41 ± 0. 1 1 . 



face. In the right panels of Figure 8 we show contours of the 
observed R' flux ratios {white curves) and la error envelops 
{colored polygons) of the various He-like ions as a function of 
7bb and W . The figure demonstrates that the observed R' flux 
ratios of N VI, O VII, Ne IX, Mg XI, and Si XIII can be pro- 
duced in a low-density plasma sitting on the white dwarf sur- 
face if the blackbody temperature T^b ~ 7, 10, 14, 18, and 30 
kK, respectively, or at higher temperatures at greater distances 
from the white dwarf. Conversely, the figure gives the allowed 
range of distances (dilution factors) for each ion for a given 
blackbody temperature. For example, for T^b = 25 kK, the vol- 
ume of plasma in which Si XIII dominates [logr(K) « 7.0] 
could be on the white dwarf surface, while that of Mg XI 
[logr(K) w 6.8] would have to be at r « 3R^i, that of Ne X 
[log r(K) K, 6.6] would have to be at r « IQR^i, and so on for 
the lower Z ions. 

2.4. Radial velocities 

In the next component of our analysis, we used two tech- 
niques to search for orbit- and spin-phase radial velocity vari- 
ations in the X-ray emission lines of AE Aqr. 

2.4.1. Composite line profile technique 

In the first technique, similar t o tha t employed by 
iHoogerwerf. Brickhouse. & Mauch j (120041) . phase-average 
as well as orbit- and spin-phase resolved composite line pro- 
files were formed by coadding the HETG event data in veloc- 
ity space V = c(A- Ao)/A relative to laboratory wavelengths 
Ao from the Interactive Guide for ATOMDB version 1.3. The 
lines used in this analysis are those labeled in Figures 5 and 
7 shortward of 20 A. For the H-like lines, the wavelengths 
are the mean of the wavelengths of the doublets weighted by 
their relative emissivities (2:1), while for the He-like lines, the 
wavelengths are for the stronger resonance lines. The result- 
ing phase-average composite line profile is shown in Figure 
9. Fit with a Gaussian, its offset v = 25 ± 26 km s"' and width 
cr = 712 ±27 km s-i. 

Applying the composite line profile technique to the orbit- 
phase resolved composite line profiles results in the radial ve- 
locities shown in Figure 10a. Assuming that, like EX Hya 
( IHoogerwerf. Brickhouse. & Mauchell2004<) . the orbit-phase 
radial velocities foHow the motion of the white dwarf, these 
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Fig. 8.— {ay(e): R' = //(i + r) flux ratio versus logWc for N VI, O VII Ne IX , Mg XI, and Si XIII for Tbb = and logTe(K) = 6.2, 6.4, 6.6, 6.8, and 7.0, 
respectively. Observed flux ratios and errors are from Table 4 of lltoh et alj 120061 (red horizontal lines) and Table 1 of this paper {blue horizontal lines); the 
inferred electron densities and en'ors (colored vertical lines) are listed in Table 1. Colored boxes delineate the Icr error envelope of the flux ratio and lognc for 
each ion. (/)-(/): Corresponding contours of the observed R' flux ratios (white curves) and Icr error envelops (colored polygons) for nc = 1 cm"'' as a function of 
Tbb and W . 



data are well fit (x^ ■ 



1 .29 /6 = 0.22) by the sine function 



v{(i)o,h) = l+K- 



1 



sin27r0t/0, 



(1) 



' </>orb-A/2 

where A = 1/8, with 7 = 30 ± 25 km s"' and = 12 ± 
36 km s"' (solid curve in Fig. lOa), although they are slightly 
better fit (xl = 1.41/7 = 0.20) with a constant 7 = 29 ± 
25 km s"'. For future reference, we note that the Icr, 2a, and 
3(7 (Ax^ = 1.0, 2.71, and 6.63) upper limits to the orbit-phase 
radial velocity semi-amplitude K = 48, 72, and 105 km s~'. 



respectively. 

Applying the composite line profile technique to the spin- 
phase resolved composite line profiles results in the radial ve- 
locities shown in Figure 10b. Unlike the orbit-phase radial 
velocities, these data are not well fit (xl = 18.3/7 = 2.6) by a 
constant, but they are well fit (xi = 2.85/5 = 0.57) by the sine 
function 



v((Aspin) 



fA/2 



sin 47r(0-(/)o) £/(/), (2) 
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Fig. 9. — MEG plus HEG composite line profile of AE Aqr (filled circles 
with error bars) and best-fitting Gaussian function {solid curve). 
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Fig. 10. — Orbit- and spin-phase radial velocities of the X-ray emission 
lines of AE Aqr. Four panels show the data (filled circles with error bars), 
best-fitting sine function (solid curve), and 7 velocity (dotted line) for (a) and 
(b) the composite line technique, (c) the cross correlation technique, and (d) 
the cross correlation technique using the boost-strapped template spectrum. 
Shaded region in the upper panel is the la error envelope of the expected 
white dwarf orbit-phase radial velocity variation. 



where A = 1 /8, with 7 = 29 ± 25 km s'^ = 149 ± 38 km 
and (po = 0.060 ± 0.021 (soUd curve in Fig. lOb). 

2.4.2. Cross correlation technique 

The composite line profile technique employed above uti- 
lizes the strongest isolated emission Unes in the HETG spec- 
trum of AE Aqr, ignoring the many weaker and often blended 
spectral features shown in Figures 5 and 7. In an attempt to 
reduce the size of the error bars on the derived spin-phase 
radial velocities, a cross correlation technique was tested. 
To accomplish this, spin-phase resolved spectra were formed 
by adding the HETG event data in wavelength space using 
bins of constant velocity width Av =100 km s~' (specifically, 
A = 3 .000, 3 .001 , 3 .002, . . .25 .002 A). In the absence of an ob- 
vious template against which to cross correlate the resulting 
spin-phase resolved spectra, the spectrum from the first spin 
phase bin, (^spin = ± A/2, was used as the template. The re- 
sulting spin-phase radial velocities, shown in Figure 10c, are 
very similar to those derived using the composite line tech- 
nique (Fig. 10/7), but the error bars are smaller by a factor 
of approximately 40%. These data are reasonably well fit 
{xl = 4.75/4 = 1.19) by equation 2 with 7 = 22 ± 15 km s"\ 
K=161± 23 km s"', and (fx,) = 0.045 ± 0.012 (solid curve in 
Fig. 10c). Note that, given the manner in which this result was 
derived, the 7 velocity is now relative that of the template, the 
spectrum from the first spin phase, for which the composite 
line profile technique gave a radial velocity v = -1 ± 57 km s"' 
(i.e., consistent with zero). 

Given this fit to the radial velocities of the spin-phase re- 
solved spectra, it is possible to produce a spin-phase aver- 
age spectrum that accounts for (removes the effect of) the 
spin-phase radial velocities. This was accomplished by mul- 
tiplying the wavelengths of the HETG event data by a factor 

1 - v(0spin)/c, where v((/)spin) = J + K sin47r((^spin-'/'o) km s"' 
with parameters that are equal to the previous best-fit values: 
7 = 22kms"\ A'= 167kms"\ and 0o = 0.045. The spin-phase 
radial velocities derived using the resulting boot-strapped 
spin-phase average spectrum as the cross correlation template 
are shown in Figure lOd. They are very similar to the radial 
velocities derived using the "vanilla" cross correlation tech- 
nique (Fig. 10c), but the error bars are smaller by a factor of 
approximately 30%. These data are now not particularly well 
fit (xl = 10.3 /5 = 2.1) by equation 2 with 7 = 15 ± 10 km s'^ 
K=163± 15 km s"', and = 0.022 ± 0.008 (solid curve in 
Fig. lOii). The deviations from the fit appear to be consistent 
with a radial velocity amplitude that is larger on the 0-0.5 
spin-phase interval and smaller on the 0.5-1 spin-phase inter- 
val. Accordingly, equation 2 was modified to allow this addi- 
tional parameter, and the data are then well fit (xi = 0.82/4 = 
0.21)with7= 14±10kms-',/:i =206±20kms-' (valid on 
^3pij, = 0-0.5), K2 = 120 ± 20 km s"' (vaUd on (/)spin = 0.5-1), 
and (f>o = 0.023 ± 0.008 (dotted curve in Fig. lOd). 

3. SUMMARY 

As summarized below, our long, uninterrupted Chandra 
HETG observation provides a wealth of details concerning 
the X-ray emission of AE Aqr: 

1. The X-ray light curve is dominated by flares that last 
between a few hundred and a few thousand seconds, pro- 
duce increases of up to 3-5 times the baseline count rate 
(Fig. 1), and are achromatic except near their peaks (Fig. 2); 
the white dwarf spin-phase X-ray light curve is achromatic 
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and sinusoidal in shape, with a relative semi-amplitude of ap- 
proximately 16% (Fig. 3); and the X-ray pulse amplitude in- 
creases linearly with the mean count rate in the middle of the 
range, but saturates at both the low and high ends of the range 
(Fig. 4). 

2. The X-ray spectrum is that of an optically thin multi- 
temperature thermal plasma (Fig. 5); the X-ray emission lines 
are broad (Fig. 9), with widths that increase with the line en- 
ergy, from cr w 1 eV (510 km s"') for O VIII to ct « 5.5 eV 
(820 km s"') for Si XIV (Fig. 6); the X-ray spectrum is rea- 
sonably well fit by a plasma model with a Gaussian emission 
measure distribution that peaks at log T(K) = 7.16, has a width 
cr = 0.48, an Fe abundance equal to 0.44 times solar, and other 
metal (primarily Ne, Mg, and Si) abundances equal to 0.76 
times solar (Fig. 7); and for a distance d = 100 pc, the total 
emission measure EM = 8.0 x lO^-' cm"-' and the 0.5-10 keV 
luminosity Lx = 1.1 x 10-" erg s"'. 

3. Based on the f/(i+r) flux ratios of the He a triplets 
of N VI, O VII, Ne IX measured by Itoh et al. in the XMM- 
Newton Reflection Grating Spectrometer spectrum, and those 
of O VII, Ne IX, Mg XI, and Si XIII in the Chandra HETG 
spectrum, the electron density of the plasma increases with 
temperature by over three orders of magnitude, from ~ 6 x 
10'" cm-3 for N VI [log 7;(K) « 6.2] to He « 1 x lO'"* cm'^ for 
Si XIII [logre(K) « 7.0] (Table 1 and Fig. %a-e), and/or the 
plasma is significantly affected by photoexcitation (Fig. 8/- 

;■)■ 

4. The radial velocity of the X-ray emission lines varies 
on the white dwarf spin phase, with two oscillations per spin 
cycle and an amplitude K « 160 km s~' (Fig. 10). 

4. DISCUSSION 

Over the years, two very different models have been pro- 
posed for the source of the X-ray emission of AE Aqr. On 
one hand, based on ROSAT and ASCA data, EracleousI i 1999b 
argued that the X-ray emission, including the flares, must 
occur close to the white dwarf, so that the gravitational po- 
tential energy can heat the X-ray emitting plasma to the ob- 
serve d temperatures. Similarly , based on Ginga and ASCA 
data, Choi. Dotani. & Agrawall (119991) argued that the X-ray 
emission, both persiste nt and flare, orig inates within the white 
dwarf magnetosphere: IChoi & Dotanil (12006) came to similar 
conclusions based on XMM-Newton Optical Monitor (OM) 
and European Photon Imaging Camera (EPIC) data. On 
the other hand, based on XMM-Newton RGS data, lltoh et alJ 
(|2006) argued that the f/(i + r) flux ratios of the He a triplets 
of N VI, O VII, and Ne IX are consistent with a plasma with an 
electron density ^ 10" cm"^ and, given the observed emis- 
sion measure, a linear scale / w (2-3) x 10"^ cm. Because this 
density is orders of magnetic less than the conventional esti- 
mate for the post-shock accretion column of a magnetic CV, 
and because this linear scale is much larger than the radius of 
the white dwarf, these authors argued that the optically thin 
X-ray-emitting plasma in AE Aqr is due not to accretion onto 
the white dwarf, but to blobs in the accretion stream, heated 
to X-ray emitting temperatu res by the propel ler action of the 
white dwarf magnetic field. Ilkhsanovl (|2006) has taken issue 
with some of the details of this model, arguing that the de- 
tected X-rays are due to either (1) a tenuous component of the 
accretion stream or (2) plasma "outside the system," heated 
by accelerated particles and/or MHD waves due to a pulsar- 
like mechanism powered by the spin-down of the magnetic 
white dwarf. The presence of non-thermal particles in AE 



Aqr is suppo rted by the o b served TeV 7-rays and the recent 
discovery by iTerada et al.l ( 120081) of a hard, possibly power- 
law, component in the Suzaku X-ray spectrum of AE Aqr 
Finally, Venter & Meintjes (2007) have proposed that the ob- 
served unpulsed X-ray emission in AE Aqr is the result of a 
very tenuous hot corona associated with the secondary star, 
which is pumped magnetohydrodynamically by the propeller 
action of the white dwarf magnetic field. As we argue below, 
the results of our Chandra HETG observation of AE Aqr — 
particularly the orbit-phase pulse time delays, the high elec- 
tron densities and/or high levels of photoexcitation implied 
by the He a triplet flux ratios, and the large widths and spin- 
phase radial velocities of the X-ray emission lines — are not 
consistent with an extended, low-density source of X-rays in 
AE Aqr, but instead support earlier models in which the dom- 
inant source of X-rays is of high density and/or in close prox- 
imity to the white dwarf. 

Consider first the systemic velocity of the X-ray emis- 
sion lines. We variously measured this quantity to be v = 
-1 ±33 km s~' from the Lyman a emission lines (§2.3.1), 
V = 25 ± 26 km s"' from the composite line profile (Fig. 9), 
7 = 29 ± 25 km s"' from the composite line profile fit to the 
radial velocities (Fig. 10b), and 7 = 22 ± 15 km s"' from 
the cross correlation fit to the radial velocities (Fig. 10c). 
Optimistically assuming that these measurements are neither 
correlated nor strongly affected by systematic effects, the 
weighted mean and standard deviation of the systemic veloc- 
ity of the X-ray emission lines 7x = 21 ± 1 1 km s"'. Relative 
to the systemic velocity of the optical emission (absorption) 
lines, 70 « -38 ± 9 (-64 ± 1 1) km s"' (Robinson et al. 199 It 
Reinsch & Beuermann 1994; Welsh et al. 1995; Casares et al] 
19961: IWatson et alJ 120061) . the X -ray emission lines are 
redshifted by Av = 7x - 70 « 59 ± 14 (85 ± 16) km s"'. 
This result is to be compared to the free-fall velocity Vff = 
(2GM„d//?wd)'''^ ~ 5500 km s"' onto the surface of the white 
dwarf, the infall velocity v'i,i < Vff/4 « 1375 km s~' below 
the putative standoff shock, and the gravitational redshift 
Av = GM„i/R„iC = 5 1^}[ km s~' from the surface of the white 
dwarf with a mass M„d = 0.8 ±0.1 M0 and radius R„i = 
(7.0 =F 0-8) X 10^ cm ( Robinson et al. 1991; Welsh et al. 19931 
Reinsch & Beuermann 1994: .Welsh et al..l995; Casares et aH 
1996t IWatson et al. 2006). Although the determination and 
interpretation of systemic velocities is fraught with uncertain- 
ties, the systemic velocity of the X-ray emission lines of AE 
Aqr is consistent with the gravitational redshift from the sur- 
face of the white dwarf, and hence with a source on or near 
the surface of the white dwarf, rather than a more extended 
region within the binary. 

Supporting the proposal that the X-ray emission of AE Aqr 
is closely associated with the white dwarf is the fact, es- 
tablished previously by Mauche (2006), that the X-ray spin 
pulse follows the motion of the white dwarf around the bi- 
nary center of mass, producing a time delay Af = 2.17 ± 0.48 
s in the arrival times of the X-ray pulses. Given this re- 
sult, it is somewhat surprising that the radial velocity of the 
X-ray emission lines does not appear to vary on the white 
dwarf orbit phase: the measured orbit-phase radial velocity 
semi-amplitude K = 12 ± 36 km s"', whereas the expected 
value K = 27rAfc/forb = 115 ± 25 km s"' (Fig. 10a). The 
difference between the measured and expected radial veloc- 
ity semi-amplitudes is 103 ±44 km s"', which differs from 
zero by 2.3cr. While this discrepancy is of some concern, we 
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showed above that the radial velocity of the X-ray emission 
lines varies on the white dwarf spin phase (Fig. lOb-d), which 
argues strongly for a source of X-rays trapped within, and ro- 
tating with, the magnetosphere of the white dwarf. 

In contrast, it seems clear that our result for the radial ve- 
locity of the X-ray emiss ion lines i s not consistent with the 
proposal, put forward by Itoh et al.l (120061) . that the accretion 
stream is the dominant source of X-rays in AE Aqr. For the 
system parameters of AE Aqr, the accretion stream makes its 
closest approach to the white dwarf at a radius r„^in « 1 x 10'" 
cm and a velocity v^ax ~ 1500 km s"' . Accounting for the bi- 
nary inclination angle / = 60°, the predicted accretion stream 
radial velocity amplitude Vn,axSin/ ~ 1300 km s"', whereas 
the 3a upper limit to the orbit-phase radial velocity semi- 
amplitude = 105 km s"'. Clearly, the accretion stream, if 
it follows a trajectory anything like t he inhomogeneous dia- 
magne tic accretion flow calculated by IWvnn. King. & Hornd 
(Il997h . cannot be the dominant source of X-rays in AE Aqr. 

An additional argument against the accretio n stream 
(lltoh et al.l 12006 ^. plasma "outside the system" (llkhsanovl 
I2006h. a ho t corona associated with the secondary star 
(I Venter & M eintjes 2007), or any other extended source of 
X-rays in AE Aqr is the high electron densities logng ~ 10.8, 
11.4, 12.3, 13.0, and 14.1 inferred from the //(/+r) flux ratios 
of the He a ti'iplets of N VI, O VII, Ne IX, Mg XI, and Si XIII, 
respectively. Given the differential emission measure distri- 
bution dEM/d\ogT « 6.6 x 10^^^ exp[-(logr-logro)V2cr^] 
with logro(K) « 7.16 and a « 0.48, the N VI, O VII, Ne IX, 
Mg XI, and Si XIII He a triplet emissivity-weighted emis- 
sion measure EM = J{dEM /dlogT)edlogT / j edlogT « 
[0.9,1.6,3.3,4.8,5.8] x lO'^^ ^m-^ and the linear scale / = 
{EM/nlfl^ « 1.4 X 10'°, 6.7 x 10^ 2.0 x 10^ 7.4 x 10^ and 
1.4 X 10^ cm or 20, 10, 3, 1, and 0.2 R„i, respectively. Al- 
though even the observed Si XIII f/(i+r) flux ratio can be 
produced in a low-density plasma suffering photoexcitation 
by an external radiation field, this requires both high black- 
body temperatures (Fbh > 25 kK) and large dilution factors 
(W < j), hence small-to-zero distances above the white dwarf 
surface. We conclude that, if the bulk of the plasma near the 
peak of the emission measure distribution is not of high den- 
sity, it must be in close proximity to the white dwarf. 

4.1. White dwarf 

To investigate the possibility that the X-rays observed in 
AE Aqr are associated with the white dwarf, we consid- 
ered a simple geomet ric model, similar to that derived by 
lEracleous et alj (119941) from the maximum entropy maps of 
the UV pulse profiles measured with the Faint Object Spec- 
trograph onboard the Hubble Space Telescope. The model 
consists of a rotating white dwarf viewed at an inclination an- 
gle / = 60° from the rotation axis z with two bright spots, cen- 
tered 20° above and below the equator and separated in lon- 
gitude by 180° (i.e., the upper and lower spots are centered at 
spherical coordinates [6',0] = [70°, 0°] and [110°, 180°], re- 
spectively). Such a model naturally produces two unequal 
flux peaks per spin cycle, with the brighter (dimmer) peak 
occurring when the upper (lower) spot is pointed toward the 
observer. Instead of the optically thick assumption applied in 
the optical and UV, we assume that the X-ray emitting spots 
are optically thin. In this case, the X-ray flux modulation 
is produced solely by occultation by the body of the white 
dwarf, so the brightness of the lower spot must be reduced by 
approximately 30% to avoid producing a second peak in the 
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Fig. 1 1 . — Spin-phase (a) light curves and (b) radial velocities for the two- 
spot model of AE Aqr. Blue, red, and black curves are for the upper spot, the 
lower spot, and the total flux, respectively. Flux is relative to the total mean 
and velocities assume I'rot = 1330 km s"' . Upper graphic shows a schematic 
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X-ray light curve. Assuming that the brightness distribution 
of the spots is given by a Gaussian function e\p{-Sd^ /2a^), 
where 69 is the polar angle from the center of the spot and the 
spot width cr = 30° (hence FWHM = 70° and the fractional 
emitting area / = J exp(-9^ /2a^) sin 9dd « 0.25, justifying 
our choice for the fiducial value of this quantity in §2.3.2), 
the light curves of the upper spot, the lower spot, and the to- 
tal flux are as shown in Figure 11a. Smaller spots produce 
squarer light curves, while larger spots produce lower relative 
oscillation amplitudes; the relative pulse semi-amplitude of 
the model shown, (/max-^min)/(^max+^iin) = 16%, is Consistent 
with observations (Fig. 3a). The flux-weighted mean radial 
velocities of the emission from the upper spot, the lower spot, 
and the total flux are as shown in Figure lib for an assumed 
white dwarf rotation velocity Vrot = 27r7?wd/^'33 = 1330 km s"' . 
As observed, the predicted radial velocity of the total flux 
goes through two oscillations per spin cycle, although the two 
peaks are not equal in strength, the red-to-blue zero velocity 
crossings do not occur near 0spin = 0.25 and 0.75, and the max- 
imum velocity amplitude is greater than observed (Fig. 10^7- 
d). Although it is possible to remedy these deficiencies by 
using two bright spots that are more nearly equal in bright- 
ness, this significantly reduces the pulse amplitude of the to- 
tal flux. Furthermore, the model predicts that the width of 
the X-ray emission lines should be narrower near 0spin = 
and 0.5 and broader near (/)spin = 0.25 and 0.75, which is not 
corroborated by the data. Finally, the model predicts that the 
radial velocity amplitude of the total flux is higher (lower) 
on the 0.75-1.25 (0.25-0.75) spin phase interval, whereas the 
boot-strapped cross correlation technique provides evidence 
that this is the case approximately one quarter of a cycle later: 
on the 0-0.5 (0.5-1) spin phase interval (ignoring the small 
phase offset). 

4.2. Accretion Columns 
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To place their re sults in a theoretical framework, 
lEracleous et alj ( 119941) interpreted the hot spots derived from 
their maximum entropy maps as the result of reprocessing of 
X-rays emitted by the post-shock gas in the accretion columns 
of a magnetic white dwarf. For an assumed point source of il- 
lumination, they found good fits to the UV pulse profiles if 
the angle between the spin axis and the magnetic axis — the 
magnetic colatitude — /3 = 76° ±1°, the peak spot tempera- 
ture 7;„ax = 26 ± 2 kK, the height of the illuminating source 
above the white dwarf surface H /R„d = 3 ± 1 , and the X-ray 
luminosity Lx = 3 x lQ^^(H/3)^(ri/Q.5)~^ erg s"', where rj is 
the efficiency of conversion of accretion luminosity into 0.1- 
4 keV X-rays (note that our fiducial value for this quantity is 
larger by an order of magnitude than that assumed by Era- 
cleous et al. because the mean temperature of the X-ray spec- 
trum of AE Aqr is uncharacteristically low). These results are 
troubling for two reasons. First, as noted by Eracleous et al., 
the shock height required to produce the size of the UV spots 
is quite large (comparable to the corotation radius). Second, 
the X-ray luminosity required to heat the 26 kK spots by re- 
processing is more than two orders of magnitude greater than 
observed in the 0.5-10 keV bandpass. Adding to these prob- 
lems, we found that such a model reproduces neither the X- 
ray pulse profile nor the radial velocity variation observed in 
AE Aqr. First, the predicted pulse profile has a square wave- 
form, with the X-ray flux rapidly doubling (halving) when one 
of the illuminating spots emerges from (is hidden behind) the 
body of the white dwarf. Second (if the illuminating spots are 
not equal in brightness), the predicted radial velocity variation 
goes through only one oscillation per spin cycle. 

Given these problems, we modeled the accretion columns 
as two uniform emission volumes contained within a polar 
angle of the magnetic axis and radii r/R„i - 1 = 0-h. For 
an inclination angle / = 60°, the pulse profiles and radial ve- 
locity variations were investigated for magnetic colatitudes 
P = 60°-90°, opening angles A0 = 15°-90°, shock heights 
h = 0-3, and infall velocities Vi„faii = O-Vff = 0-5500 km s"', 
with and without intensity weighting by a Gaussian func- 
tion e\p(-S9- /2<T^). Although the parameter space is huge, it 
did not seem possible to reproduce the observed X-ray pulse 
profile, let alone the radial velocity variation, with such a 
model, unless it looked very much like the two-spot model 
discussed in the previous section (i.e., f3 « 70°, a w 30°, 
h « 0, and vinfaii < 0.01 Vff): even modest infall velocities pro- 
duced model radial velocity variations with systematic veloc- 
ities that are greater (in absolute magnitude) than observed, 
while large shock heights and/or opening angles and/or Gaus- 
sian widths produced small relative pulse amplitudes, and, 
for much of the parameter space, flux maximum occurs at 
spin phases ip^pin = 0.25 and 0.75, when both of the accretion 
columns are visible on the limb of the white dwarf. 

4.3. Magnetosphere 

We next considered the possibility, argued by 
IChoi. Dotani. & Agrawall (Il999h . that the X-ray emis- 
sion of AE Aqr, both persistent and flare, is due to plasma 
associated with the magnetosphere of the white dwarf. We 
assume that the magnetosphere is filled with plasma stripped 
off of the accretion stream, which, as noted above, makes its 
closest approach to the white dwarf at a radius r^i^ w 1 x 10'" 
cm and a velocity v^ax ~ 1500 km s"' . The magnetic field of 
the white dwarf can control the motion of an ionized com- 
ponent of the stream if the magnetic pressure B^/Stt [where 



B = B*(/?wd/'"min)^ for a dipole field with a surface magnetic 
field strength B^,] is greater than the ram pressure pv^ax' hence 
if > (87r^mHne)'/'v™ax(r,n,„/-Rwd)^ « 7(«e/10i^^ cm-3)i/2 
MG. If the kinetic energy of the accretion stream is 
thermalized in a strong shock at its closest approach 
to the white dwarf, it will be heated to a temperature 
Tmax = 3/imHV^ax/16A: « 3.3 X 10^ K or approximately twice 
the peak temperature of the emission measure distribution. 
Hence, the magnetosphere plausibly could be the source of 
X-rays in AE Aqr, independent of any other source of energy, 
ifit is fed at a rate M = 2Lx/vLx ~ 9.8 x lO'^^gs"'. 

The magnetosphere is an attractive source of X-rays in AE 
Aqr in as much as it would naturally supply a range of densi- 
ties, linear scales, and velocities, as required by the data. On 
the other hand, we found that the widths of the X-ray emis- 
sion lines increase from w 1 eV (510 km s~') for O VIII 
to cr « 5.5 eV (820 km s"') for Si XIV (Fig. 6), whereas 
the line widths predicted by this model are nominally much 
larger: for plasma trapped on, and forced to rotate with, the 
white dwarf magnetic field, the projected rotation velocity 
v,-otsin/ = 27rrsin//f33, which varies from 1150 kms"' forr = 
^wd = 7x10** cm (as in the two-spot model) to 16,500 km s~' 
for r = r^iin = 1x10'" cm. However, the observed value for the 
lines widths should be near the lower end of this range, since, 
for a magnetosphere uniformly filled with plasma with an in- 
ward (radial) velocity Vi- oc Vff = (2GMwd/'")'''^, the plasma 
density p oc r~^l'^ and the X-ray emissivity e oc oc r"^. 

5. CONCLUSION 

Of the simple models considered above for the source of 
the X-ray emission of AE Aqr, the white dwarf and mag- 
netosphere models are the most promising, while the ac- 
cretion column model appears to be untenable: it fails to 
reproduce the Chandra HETG X-ray light curves and ra- 
dial velocities, and it requires an X-ray luminosity that is 
more than two orders of magnitude greater than observed 
in the 0.5-10 keV bandpass to heat the UV hot spots by 
reprocessing. A more intimate association between the X- 
ray and UV emission regions could resolve this problem, 
and it is interesting and perhaps important to note that such 
a situation is naturally produced by a bombardment and/or 
blobby accretion solution to the accretion flow. In par- 
ticular, if the cyclotron-balanced bombardment solution ap- 
plies to AE Aqr, it would naturally produce comparable X- 
ray and UV spot sizes and luminosities, comparable rela- 
tive pulse amplitudes, comparable orbit-phase pulse time de- 
lays, the observ ed tight correl ation between the X-ray and 
UV light curves (lMauchell2009h . and a n accretion region that 
is heated to the observe d T ^ 10^ K (Woel k & BeuermannI 
[1992, 1993, 1996; Fische r & Beuer mann 20Qil). Such a so- 
lution to the accretion flow applies if (1) the specific accre- 
tion rate (;« = M /A) onto the white dwarf is sufficiently low 
that the accreting plasma does not pass through a hydrody- 
namic shock, but instead is stopped by Coulomb interactions 
in the white dwarf atmosphere, and (2) the magnetic field of 
the white dwarf is sufficiently strong that the accretion lumi- 
nosity can be radiated away by a combination of optically 
thin bremsstrahlung and line radiation in the X-ray wave- 
band and optically thick cyclotron ra diation in the infrared 
and op tical wavebands. According to iFischer & BeuermannI 
( 120011 equating from eqn. 17 with that from eqn. 20), the 
limit for the validity of cyclotron-balanced bombardment is 



12 



Mauche 



ni < 1.19;!;o 47 ^ Bij-^ g cm ^ s ' for a white dwarf mass 
M„d = 0.8 ±O.2M0 and magnetic field strength B-i = B*/10^ 
G. 

Are these conditions satisfied in AE Aqr? Perhaps. First, 
assuming that the observed X-ray luminosity is driven by 
accretion onto the white dwarf, the accretion rate M = 
i-x^wd/GMwd « 7.3 X lO'-' g s"', hence the specific accre- 
tion rate m « 4.7 x 10"^ (//0.25)"' g cm"^ s"'. Second, the 
strength of the magnetic field of the white dwarf in AE Aqr 
is uncertain, but based on the typical magnetic moments of 
intermediate polars (10^^ G cm-') and polars (10-'"' G cm-'), it 
should lie in the range 0.3 MG < < 30 MG. More spe- 
cific estimates of this quantity include < 2 MG, based on 
evolutionary considerations (Meintjes 2002); B^, ^ 1-5 MG, 
based on the low levels of, and upper limi t s on, the circu- 
lar polariz ation of AE Aqr (|Croppeii[T986t IStockman et al.l 
[1992; Besk rovnaya et al although these values are very 
likely lower limits, since these authors do not account explic- 
itly for AE Aqr's low accretion luminosity and bright sec- 
ondary); and Bi, ^ 50 MG, based on the (probably incor- 
rect) assumption that the s pin-down of the white dwarf is due 
to the pulsar mechanism (IIkhsanovlll998l) . The above limit 
for the validity of cyclotron-balanced bombardment requires 
B* > 7![;^(//0.25)"°-^^ MG, which is probably not inconsis- 
tent with the circular polarization measurements and upper 
limits. Conversely, the validity of the bombardment solution 
in AE Aqr could be tested with a more secure measurement 
of, or lower upper limit on, the magnetic field strength of its 
white dwarf. 

Where do we go from here? Additional constraints on 
the global model of AE Aqr, and in particular on the con- 
nection between the various emission regions, can be sup- 
plied by simultaneous multiwavelength observations, includ- 
ing those obtained during our 2005 multiwavelen gth cam- 
paign, the analysis of which is at an early stage ("Mauche* 
12009). TeV observations with the Whipple Obser vatory 10- 
m tel escope obtained over the epoch 1991-1995 CLang et al.i 
119981) and the Major Atmospheric Gamma-Ray Imaging 
Cherenkov (MAGIC) 17-m telescope obtained during our 
campaign (ISidro et al.l 120 08") provide only (low) upper lim- 
its on the TeV 7-ray flux from AE Aqr, so it would be use- 
ful to unambiguously validate or invalidate earlier reports of 
detections of TeV flux from this source; in particular, a ded- 
icated High Energy Stereoscopic System (HESS) campaign 
of observations, preferably in conjunction with simultaneous 
X-ray observations, is badly needed. Similarly, it would be 
useful to independently validate or invalidate the existence of 
the high-energy emission detected by Terada et al. (2008) in 
the Suzaku X-ray spectrum of AE Aqr, and to establish un- 



ambiguously if it is thermal or nonthermal in nature. Going 
out on a limb, it is our bet that AE Aqr is not a TeV 7-ray 
source and that the high-energy X-ray excess is thermal in na- 
ture, specifically that it is due to the stand-off shock present 
when the specific mass accretion rate onto the white dwarf is 
large (e.g., during flares) and the bombardment solution is no 
longer valid. Additional high resolution X-ray spectroscopic 
observations are also warranted, particularly to confirm or re- 
fute the high electron densities inferred from the flux ratios 
of the Ne IX, Mg XI, and Si XIII He a ti'iplets; it is, after all, 
a priori unlikely that each of these ratios lies on the knee of 
the theoretical f /(i + r) flux ratio curves, and the trend of in- 
creasing density with increasing temperature is opposite to 
what is expected from a cooling flow (although it is consis- 
tent with the opposite, e.g., adiabatic expansion). It would be 
extremely h elpful to apply the various Fe L-shell de nsity di- 
agnostics (Mauc he. Liedahl. & Foumierf200lLl2003ll20"05D to 
AE Aqr, since they are typically less sensitive to photoexci- 
tation, but the Fe XVII 17.10/17.05 line ratio for one is com- 
promised by blending due to the the large line widths, and the 
other line ratios all require higher signal-to-noise ratio spec- 
tra than currently exists. It would be extremely useful to (1) 
determine the radial velocities of individual X-ray emission 
lines, rather than all the lines together; (2) resolve the indi- 
vidual lines on the white dwarf spin and orbit phases, (3) in- 
vestigate the flare and quiescent spectra separately, (4) resolve 
individual flares, and (5) include the Fe K lines into this type 
of analysis. Such detailed investigations are beyond the ca- 
pabilities of Chandra or any other existing X-ray facility, but 
would be possible with future facilities such as IXO that sup- 
ply an order-of-magnitude or more increase in effective area 
and spectral resolution. 
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